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Genetic basis of phenotypic differences in individuals is an important area in biology and personalized medicine. 
Analysis of divergent Saccharomyces cerevisiae strains grown under different conditions revealed extensive 
variation in response to both drugs (e.g., 4-nitroquinoline 1 -oxide [4NQO]) and different carbon sources. 
Differences in 4NQO resistance were due to amino acid variation in the transcription factor Yrrl. Yrrl YJM789 
conferred 4NQO resistance but caused slower growth on glycerol, and vice versa with Yrrl S96 , indicating that 
alleles of Yrrl confer distinct phenotypes. The binding targets of Yrrl alleles from diverse yeast strains varied 
considerably among different strains grown under the same conditions as well as for the same strain under 
different conditions, indicating that distinct molecular programs are conferred by the different Yrrl alleles. Our 
results demonstrate that genetic variations in one important control gene (YRR1), lead to distinct regulatory 
programs and phenotypes in individuals. We term these polymorphic control genes "master variators." 
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How genetic variation produces different phenotypes and 
the ability to adapt to distinct environmental challenges 
is an important question in biology. Polymorphisms in 
coding regions and regulatory elements have been shown 
to associate with changes in individual morphological 
phenotypes of multicellular organisms and survival char- 
acteristics in microorganisms (Rebeiz et al. 2009; Chan 
et al. 2010; Gerke et al. 2010; Davidson 2011). However, 
the extent to which variations in a single locus can generate 
entirely distinct regulatory programs and phenotypes is 
not well understood. 

Alterations in gene regulatory components have been 
shown to vary in a range of organisms and can lead to 
distinct phenotypes. Loss of activity of transcription 
factors (TFs) can lead to phenotypic differences,- for 
example, the loss of Pitxl expression in stickleback fish 
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leads to different body morphologies (Chan et al. 2010; 
Infante et al. 2013), while variation in some TFs (e.g., 
Imel, Rmel, and Rsfl) can alter sporulation efficiency in 
different Saccharomyces cerevisiae strains (Gerke et al. 
2009). Additionally, differences in gene expression and 
TF-binding sites have been shown to vary considerably 
between closely related species and individuals, although 
how these differences affect phenotypes is often un- 
known (Borneman et al. 2006; Zheng et al. 2010). Al- 
though we know that phenotypic changes can result from 
total loss of activity of TFs, it remains unclear whether 
allelic variations in a single TF can mediate broad dif- 
ferences in molecular responses and complex phenotypes 
through network rewiring. Such changes would be 
expected to have important implications in the rapid 
evolution of complex phenotypes and in allowing the 
adaptation of organisms to new environments. 

© 2014 Gallagher et al. This article, published in Genes & Develop- 
ment, is available under a Creative Commons License (Attribution- 
NonCommercial 3.0 Imported), as described at http://creativecommons. 
org/licenses/by-nc/3.0/. 
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To better understand how phenotypes vary and evolve 
among different individuals within a population or among 
populations within a species, we analyzed different yeast 
strains to find proteins that can regulate the phenotypic 
response to the environment. Natural yeast strains are 
more genetically diverse than humans (Wei et al. 2007) 
and have been shown to vary considerably from labora- 
tory strains in their response to different environmental 
conditions such as growth on distinct carbon sources and 
incubation in the presence of compounds and drugs (Liti 
et al. 2009). The loci mediating some phenotypic differ- 
ences among different yeasts have been mapped (Gerke 
et al. 2009; Ehrenreich et al. 2010; Zheng et al. 2010), 
although the number of genes that have actually been 
identified is very limited. 

We mapped the genetic basis of resistance/sensitivity 
to 4-nitroquinoline 1 -oxide (4NQO), a bulky alkylating 
agent, in two yeast strains with very different responses 
to this drug and identified a gene, YRR1, whose variation 
strongly contributed to this phenotypic difference. En- 
hanced fitness in one condition resulted in diminished 
fitness in another. When cells were able to grow well in 
the presence of 4NQO, they then grew poorly in glycerol. 
The phenotypic difference was largely associated with 
changes in a single amino acid in Yrrl that was a potential 
site of phosphorylation. Importantly, Yrrl bound differ- 
ent gene targets in different strain backgrounds, indicat- 
ing that variation in a single TF can rewire transcriptional 
networks. Thus, natural polymorphisms in a single reg- 
ulatory gene can cause broad molecular changes that 
result in distinct phenotypes in individuals, a phenome- 
non that has not been directly observed previously. To 
distinguish this kind of protein sequence variation from 
deleterious mutations, we term such genes whose natural 
polymorphisms that cause broad molecular changes mas- 
ter variators. 

Results 

Different growth phenotypes among different 
yeast strains 

To explore the role of natural genetic variation in causing 
phenotypic differences, we first analyzed the growth of 
five haploid S. cerevisiae strains from diverse genetic and 
ecological backgrounds under a variety of environmental 
conditions. The strains examined were a common labo- 
ratory strain, S96 (derived from S288c background); two 
clinical isolates, YJM789 and YJM339 (McCusker et al. 
1994); and two vineyard isolates, AWRI1631 (Borneman 
et al. 2008) and RM11 (Mortimer et al. 1994; Brem et al. 
2002). The phenotypes analyzed included relative growth 
in different carbon sources (glucose, ethanol, galactose, 
and glycerol) and growth in the presence of different drugs 
(4NQO, fluconazole, hydroxyurea, and nystatin) (Fig. 1A). 
Fluconazole and nystatin both affect sterol metabolism, 
and YJM789 sensitivity to these chemicals has been 
recently linked to polymorphisms in the multidrug trans- 
porter Pdr5 (Guan et al. 2010). The genetic variation of 
response to hydroxyurea remains unknown. The five strains 
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Figure 1. Extensive variation in yeast growth under different 
conditions and mapping of 4NQO response. [A) Tenfold serial 
dilution of various strains of yeast spotted onto plates contain- 
ing YPD, YPD + 0.2 [ig/mL 4NQO, YP + 3% glycerol, YP + 2% 
galactose, 2% ethanol, YPD + 5 |xg/mL fluconazole, YPD + 200 mM 
hydroxyurea, and YPD + 2 fxg/mL nystatin and grown for 1-3 d at 
28°C before being photographed. (B) Tenfold serial dilution of 
S96, YJM789, and six segregants spotted onto plates containing 
YPD and 0.2 |xg/mL 4NQO and grown for 1-3 d at 28°C before 
being photographed. (C) Genome-wide analysis study identified 
a region on chromosome 15 with a LOD score of 25.50 for the 
95% Bayes credible interval 639,300-640,218. 



exhibited distinct growth characteristics under each of 
the different conditions (Fig. 1A), and each strain had a 
unique pattern of phenotypes. For example, relative to 
other yeasts, RM11 grew well under fluconazole, and 
YJM789 strains grew well in the presence of 4NQO; S96 
cells grew slightly better on glycerol relative to YJM789 
cells. Because the 4NQO resistance/sensitivity pheno- 
type showed the most variation between the strains, we 
chose to map the genetic basis of this phenotype using 
the YJM789 and S96 strains, which exhibited the most 
extreme phenotypes. 

To map the loci responsible for the 4NQO response 
variation, we used two previously sequenced yeast strains 
(YJM789 and S96) and a large collection of 125 haploid 
segregants that resulted from their cross (Wei et al. 2007; 
Mancera et al. 2008); the genetic content of each of these 
segregants was mapped at high resolution with high- 
density microarrays. As shown in Figure IB and Supple- 
mental Figure SI, variation in the 4NQO response among 
the different segregants was observed. The response of 
segregants was scored and used as a quantitative trait for 
linkage analysis between the individual trait and genetic 
markers (see the Supplemental Material; Zheng et al. 
2010). A strong peak associated with 4NQO response was 
mapped to the right arm of chromosome 15 with a peak 
logarithm (base 10) of odds (LOD) score of 25.5 (Fig. 1C). 
Four other identified genomic locations, on chromosomes 
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1, 4, 5, and 7, had LOD scores outside the 95% confidence 
interval (Supplemental Table 1) and thus were not in- 
vestigated further. 

Yrrl mediates distinct responses to 4NQO 
and growth on glycerol in different yeasts 

Because of the large number of segregants analyzed, the 
region associated with 4NQO response could be mapped 
to a very narrow interval (Fig. 2A). Within this region was 
YRR1, a gene encoding a TF. A strain containing a hyper- 
morphic allele of YRR1 was previously found to be 
resistant to 4NQO (Cui et al. 1998), but the variation of 
YRR1 among different yeasts and its potential role in 
conferring differential 4NQO resistance had not been 
studied previously. To determine whether YRR1 was re- 
sponsible for mediating variation in the 4NQO response, 
we swapped the YRR1 alleles (termed YRR1 YJM789 and 
YRR1 S96 ) between the two genetic backgrounds. YRR1 S96 
and YRRl 1 ^™ 789 alleles on centromeric plasmids were 
transformed into both wild-type and yrrlA mutant 
strains. Importantly, when either the S96 wild- type strain 
or its yrrl A derivative was expressing YRRl 1 ^™ 789 , it 
showed a strong increase in 4NQO resistance compared 
with S96 carrying an empty plasmid (Fig. 2B, rows 2, 7). 
As expected, introduction of YRR1 Y]M789 into YJM789 
yrrl A cells restored 4NQO resistance (Fig. 2B, rows 1, 
4 vs. control S96 rows 2, 8; Supplemental Fig. S2). Wild- 
type YJM789 cells containing the YRR1^ M789 allele were 
4NQO-resistant, also as expected. YJM789 yrrlA strains 
containing the YRR1 S96 allele as well as cells lacking any 
plasmid had intermediate sensitivity to 4NQO, indicat- 
ing that YRR1 is the major contributor to 4NQO re- 
sistance in this strain (Fig. 2A,B). We can also conclude 
that the YRRl^™ 789 allele is dominant over the YRR1 S96 
allele. The strong differences in growth observed on plates 
were also evident using quantitative growth assays of cells 
grown in liquid cultures containing 4NQO (Fig. 2C). 

In our initial screens, we found that S96 cells grew 
better on glycerol medium than YJM789 cells. We there- 
fore examined whether YRR 1 was involved in this pheno- 
typic difference. Remarkably, we found that YJM789 
strains (both wild type and yrrl A) exhibited a modest, 
but reproducible, increased growth on plates containing 
glycerol medium with expression of the YRR I s96 allele. 
Interestingly, S96 cells containing the YRR1 YJM789 allele 
exhibited slightly slower growth on glycerol plates. The 
results found on the plates were also reflected in the 
quantitative cell culture assays (Fig. 2C) ; YJM789 cells 
with the YRR I s96 allele grew 16% faster on glycerol 
than wild-type YJM789 cells, while the S96 cells with 
yrr1 yjm789 had a 17 o /o s i ower growth rate on glycerol 

than wild- type S96 cells. Overall, these results indicate 
that the YRR1 YJM789 allele provided increased resistance 
to 4NQO in the S96 background but reduced growth 
on glycerol in either the YJM789 or S96 background. In 
contrast, the YRR I s96 allele provided little to no resis- 
tance to 4NQO but led to increased growth on glycerol 
in either strain background. These results indicate that 
different YRR1 alleles confer distinct phenotypes and 
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Figure 2. Different alleles of YRR1 are responsible for 4NQO 
resistance and rate of growth on glycerol. (A) Ninety-five 
percent Bayes credible interval for the QTL region on chromo- 
some 15 (coordinates 636,762-645,935) of the 4NQO data. The 
LOD curve was calculated by standard interval mapping using 
EM algorithm, and the Y-axis is 10 A LOD rescaled to make the 
area under the curve within 95% interval (shaded area) equal 1. 
Associated marker names are also shown. Note that two markers 
fall on the same map position, although their physical locations 
are different (chromosome 15: 638,137-638,188, and chrl5: 
638,500-638,734). Genes are drawn to physical scale, and the 
interval for each genetic marker is noted by lines between the 
physical map in kilobases and map position in centimorgans. 
Within this interval are the following genes: PNS1 (coordinates: 
638,559-636,940) and YRR1 (coordinates: 641,993-639,561). [B] 
Tenfold serial dilution of wild-type strains YJM789 and S96 and 
yrrlA mutants carrying an empty plasmid (pGS35, a KAN R 
plasmid), pY^JRi YIM789 , or pYRRl S96 encoding the respective 
alleles of YRR1 expressed from the endogenous promoters. 
Strains were growth on YPD, YPD + 0.25 |xg/mL 4NQO, and 
glycerol with 200 |xg/mL G418 and grown for 1-3 d at 28°C before 
being photographed. (C) For cells grown in culture, the percentage 
change in doubling time for YJM789 yrrlA and S96 ynlA carrying 
an empty plasmid or the exchanged allele as compared with 
YJM789 yrrlA and S96 yrrlA carrying their respective alleles. 

that the YRR1 variations that improve growth in the 
presence of 4NQO cause reduced growth in glycerol. 

Differences in coding sequences confer phenotypic 
differences 

Inspection of the predicted protein sequences of Yrrl 
reveals that they contain multiple amino acid polymor- 
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phisms among the YJM789, S96, AWRI1631, RM11, and 
YJM339 strains. Four out of the five polymorphisms 
within the coding region were nonconservative changes 
and were located near the C terminus, and all were 
outside the putative Zn finger DNA-binding, coiled-coil, 
and activation domains (Fig. 3A). There were also four 
polymorphisms seen in the promoter regions of YRR1 YJM789 
and YRR1 S96 . Two polymorphisms downstream from the 
YRR1 gene are not transcribed (Xu et al. 2009) and were 
not further studied. To determine whether the different 
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Figure 3. Yrrl contains polymorphisms in coding regions that 
are responsible for most of the variation in 4NQO response. (A) 
Sequence alignment of Yrrl from five yeast strains. Polymor- 
phisms are indicated with Zn finger, coiled-coil, and activation 
domain boxed. [B) Immunoblot analysis of Yrrl-Myc immuno- 
precipitated from lysates from YJM789, S96, AWRI1631, RM11, 
and YJM339 strains grown in YPD, 4NQO for 2.5 h, and glycerol 
for 8.5 h. (C) The levels of Yrrl immunoprecipitated via the 
13xMyc tag from yeast were quantitated relative to Yrrl from 
YJM789 grown in YPD. An average of two to three experiments 
is shown with standard deviations. [D] Tenfold serial dilution of 
S96 yrrl A and YJM789 yrrl A carrying YRR1 with promoters 
exchanged on plasmid pGS35; strains were incubated on YPD 
and YPD + 0.25 |xg/mL 4NQO plates for 1-2 d at 28°C before 
being photographed. (£) Immunoblot of N-terminally Myc tagged 
Yrrl expressed from a plasmid with allelic promoter or trans- 
allelic promoter and immunoprecipitated from cells grown in 
YPD. Cellular extract from cells carrying Yrrl YJM789 -Myc tag was 
incubated with protein A-conjugated agarose beads to control for 
nonspecific interactions and run in lane 1, labeled Mock. 



phenotypes among strains were due to differences in the 
YRR1 promoter or in their coding sequences, several tests 
were performed. 

We first examined Yrrl protein levels in all five back- 
grounds using immunoblot analyses (Fig. 3B ; see also 
the legend). Yrrl was evident in all five strains (Fig. 3B), 
although the levels varied modestly between strain and 
condition. For cells grown on YPD and other conditions, 
the levels of Yrrl from all strains were compared with 
Yrrl YjM789 grown in YPD. The level of Yrrl S96 protein 
was 0.5 ± 0.2-fold (YPD) or less (4NQO and glycerol) than 
that of Yrrl YIM789 protein (average of three experiments) 
(Fig. 3C). There was also a similar lower level of S96 YRR1 
mRNA (0.60) compared with that of YJM789 using DNA 
microarray analyses (data not shown). However, the level 
of Yrrl did not correlate to the rate of growth in any con- 
dition or strain. Therefore, to directly determine whether 
differences in expression of YRR1 itself or coding poly- 
morphisms alter the 4NQO response, we exchanged the 
promoter regions of YRR1 from YJM789 and S96 back- 
grounds. The promoters from the two alleles of YRR1 
were cloned in front of the other YRR1 allele on a plasmid 
and expressed in both YJM789 yrrl A and S96 yrrl A. In 
S96 cells, we detected only a very small increase of 4NQO 
resistance when YRR1 S96 expression was driven by the 
yrr1 yjm789 promoter ( Fig> 3D ) Expression of the 
YRR1 Y 1 M789 allele under the control of the YRR1 S96 pro- 
moter conferred strong 4NQO resistance and was similar 
to YRR1 YJM789 under its endogenous promoter (Fig. 3D). 
Protein levels were also assessed by immunoblotting of 
Myc-tagged Yrrl proteins immunoprecipitated from 
yeast with both YRR1 alleles driven from each's endoge- 
nous promoter or the promoter from the other allele. 
There was only a slight increase of Yrrl S96 protein driven 
from the YRR1 YJM789 promoter and a slight decrease of 
Yrrl Y ™ 789 protein driven from the YRR1 S96 promoter 
(Fig. 3E). Therefore, although YRR1 expression levels may 
contribute slightly to the different phenotypes, the poly- 
morphisms in the coding region appear to be primarily 
responsible for alterations in Yrrl. 

Specific variable amino acids in Yrrl regulate 
4NQO response 

Because all of the variation among the alleles of YRR1 
was outside the DNA-binding domain, the variation in 
Yrrl function appeared to due to differences that lay in 
other domains of the protein. Longer exposures of the 
immunoblots revealed higher-molecular-weight proteins 
that were detected with anti-Myc (Fig. 3E). Mock immu- 
noprecipitations of Myc-tagged Yrrl showed no nonspe- 
cific interactions (Fig. 3E, lane 1). To investigate the 
nature of the slower-migrating bands of Yrrl, we treated 
lysates with phosphatase (calf intestinal phosphatase 
[CIP]) and immunoprecipitated Yrrl. The band with the 
same migrating pattern as the predicted Yrrl-Myc protein 
showed increased mobility, and some of the slower bands 
also shifted down to the faster-migrating band. The fastest- 
migrating band of Yrrl also showed a modest increase 
when immunoprecipitated Yrrl YIM789 was treated with 
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CIP (Fig. 4A). At least one of the post-translational 
modifications of Yrrl was sensitive to dephosphorylation 
treatment. 

We compared the five coding polymorphisms between 
Yrrl YIM789 and Yrrl S96 to examine whether any had an 
effect on Yrrl function, potentially through phosphory- 
lation. T775 and S759 were the only two variable residues 
that could be phosphorylated in Yrrl YIM789 and not in 
Yrrl S96 . To determine whether T775 has a role in Yrrl 
4NQO resistance, we made site-directed mutations in 
Yrrl YjM789 and Yrrl S96 (Fig. 4B). Serine and threonine are 
amino acids with polar uncharged side chains. When 
threonine and serine are phosphorylated, they become 
strongly negatively charged. Mutating threonine and 
serine to glutamate mimics that strong negative charge 
of phosphorylation. Mutating these residues to an ala- 
nine, a nonpolar amino acid, would prevent a negative 
charge. We tested the effect of site-directed mutation by 
mutating variable residues on plasmid-encoded Yrrl to 
glutamate (E) or alanine (A) and testing the growth of 
transformed S96 yrrlA. Mutation of amino acid T775 to 
E in both the Yrrl YJM789 and Yrrl S96 alleles increased 
4NQO resistance to even greater than the wild-type 
Yrrl YJM789 (Fig. 4B). Expression of the yrrl YJM789 T775E 
allele conferred greater 4NQO resistance than yrrl YJ]VI789 
T775A or Yrrl YJM789 . Therefore, yrrl YJM789 T775E and 
yrrl S96 T775E were hypomorphic alleles. yrrl YJM789 
T775A and yrrl S96 I775A did not change 4NQO resis- 
tance compared with the respective wild- type Yrrl al- 
leles. The yrrl Y ™ 789 S759E mutation did not signifi- 
cantly change the 4NQO resistance of the strain, while 
the yrrl YJM789 S759A resulted in greater 4NQO sensitiv- 



ity. Mutating 775 to a negatively charged amino acid 
conferred greater 4NQO resistance than the resistant 
Yrrl Y ™ 789 allele alone. We postulate that while T775 
could be the major residue of interaction, the increased 
sensitivity of yeast expressing yrrl YJA/I789 S759A suggests 
that both residues are sites of regulation. We made 
a double mutant (yrrl YJM789 S759AT775A), expressed it 
in S96 yrrlA yeast, and found that the yeast were more 
sensitive than with Yrrl^ 1 ^ 789 bearing either single muta- 
tion but not as sensitive as Yrrl S96 -expressing yeast (Fig. 4C). 

We examined two variable residues of Yrrl that are 
potential phosphorylation sites and mutated them to 
glutamate to change the charge. Only one other site in 
Yrrl varies in a way that changes the charge of the pro- 
tein. Amino acid 673 is a glutamate in Yrrl YIM789 and 
a glycine in Yrrl S96 , giving rise to different charges at this 
site. To investigate the impact of amino acid 673 between 
the different alleles of Yrrl, we mutated that site and 
tested the 4NQO resistance. S96 yrrlA yeast expressing 
yrrl YjM789 E67 3Q were as 4NQO-sensitive as yeast 
expressing Yrrl S96 (Fig. 4B). However, the converse was 
not true: Yeast expressing yrrl S96 G673E were as 4NQO- 
sensitive as Yrrl S96 -expressing yeast. We further investi- 
gated these mutants by introducing mutations at 775 
(Fig. 4D). Both the yrrl YJM789 T775E and yrrl S96 I775E 
alleles conferred strong 4NQO resistance to yeast. 
yrrl YjM789 E673G converted this allele from YJM789 to 
a 4QNO-sensitive one, and combining this mutation 
with T775E restored 4NQO resistance. The double mu- 
tants of both alleles (yrrl YJM789 E673G T775E and yrr S96 
G673E I775E) were converted to 4NQO-resistant alleles. 
Consistent with multiple post-translational modifica- 




Figure 4. Variation in a single amino acid of Yrrl Y ^ M789 contributed to 4NQO resistance. [A) Immunoprecipitation of pMyc-Yrrl Y * M789 
from S96 lysate treated overnight with CIP. [B) Tenfold serial dilution of S96 yrrl A carrying either plasmid with the indicated single 
amino acid mutant alleles of Yrrl YJM789 or Yrrl S96 on YPD and YPD + 0.3 |xg/mL 4NQO plates for 2 d at 28°C before being 
photographed. The mutations 775E, 775 A, 756E, and 756A were generated in both alleles of Yrrl. The amino acid at position 673 was 
polymorphic between the two alleles and was switched in Yrrl alleles to generate yrrl YJM798 E673G and yrrl S96 G673E alleles. S96 
yrrlA with an empty plasmid is noted with a dash at the top of each plate. (C) A double mutant of yrrl Y J M789 S796A and T775A was also 
transformed into S96 yrrl A yeast, and the resistance to 0.3 |xg/mL 4NQO was tested. (D) Double amino acid substitutions at positions 
673 and 775 as indicated were generated for each Yrrl allele, and the growth on 0.25 |xg/mL 4NQO was tested. 
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tions of Yrrl, there were multiple variable amino acids 
that regulate Yrrl's response to 4NQO. 

Assessing the effect of yul S96 I775E 
on the genomic localization 

Mutation of the variable amino acid 775 to a glutamate 
increased the 4NQO resistance above the 4NQO resis- 
tance allele Yrrl YIM789 . This mutation is in the C-termi- 
nal region of the protein and is not predicted to affect the 
Zn finger-binding domain (Fig. 3A). To test whether the 
yrrl S96 1775E allele would change the genomic location of 
the TF, we mapped the chromatin-binding sites for the 
different Yrrl alleles using chromatin immunoprecipita- 
tion (ChIP) followed by Illumina sequencing (ChlP-seq), 
and the reads were mapped back to the reference genome 
using CisGenome (Supplemental Table 2; Johnson et al. 
2007; Robertson et al. 2007; Ji et al. 2008). N-terminal 
13xMyc-tagged Yrrl alleles were expressed from plas- 
mids that were transformed into S96 yrrl A; this design 
avoided complications due to differences in genetic 
backgrounds and also ensured that the only Yrrl allele 
expressed was the one present on the plasmid. S96 yrrlA 
strains carrying either Yrrl YIM789 , Yrrl S96 , or yrrl S96 
I775E were grown in YPD only or in YPD and then shifted 
to 4NQO for 2.5 h. Sequenced reads were obtained from 
three independent biological replicate experiments, along 
with "input" controls (which gave little background, as 
described previously,- Lefrancois et al. 2009). Across the 
three strains and two conditions, a total of 329 peaks were 
identified using CisGenome (Ji et al. 2008) as potential 



binding sites of Yrrl (Supplemental Table 3). Cluster 
analysis was performed on the ChlP-seq peaks to assess 
any similarities in patterns of TF -binding. The yrrl S96 
I775E peaks from YPD-and-4NQO-treated cells clustered 
with Yrrl YJM789 peaks from 4NQO-treated cells but not 
with Yrrl S96 peaks from YPD- or 4NQO-grown cells 
(Fig. 5A). While these patterns of Yrrl binding were the 
most similar, there were additional loci bound by Yrrl S96 
I775E in 4NQO-treated cells not seen in Yrrl S96 I775E 
from YPD-grown cells. There is a group of loci bound by 
Yrrl YIM789 , but there was little to no binding by Yrrl S96 
I775E when cells were grown in YPD or 4NQO. This 
suggests that although amino acid 775 played the major 
role in regulating 4NQO responses, some other variable 
amino acids also contributed to the binding variation of 
Yrrl in response to 4NQO. 

To assess the impact of Yrrl S96 on individual loci, three 
loci were closely examined for each condition and allele. 
We saw that the mutation I775E converted Yrrl S96 to 
mimicking the behavior of Yrrl YJM789 allele in cells grown 
in 4NQO (Fig. 5B ; Supplemental Table 4). Yrrl YJM789 from 
cells grown in 4NQO (but not in YPD alone) as well as 
Yrrl s96 I775£ from cells grown in ^^er 4NQO or YPD 

alone bound upstream of YAR068w. The peaks of 
Yrrl YjM789 and Yrrl S96 I775E grown in 4NQO were simi- 
lar, while the peak from Yrrl S96 I775E grown in YPD was 
narrower. Expression of YAR068w is induced in petite cells 
(respiration-deficient yeast) (Epstein et al. 2001). The peak 
over coding region of FMP52, a gene encoding a mitochon- 
drial protein that is induced in DNA damage (Dardalhon 
et al. 2007), also shifted between alleles and conditions. 




Figure 5. Genomic locations of yrrl S96 1775 compared with Yrrl S96 and Yrrl YJM789 grown in 4NQO. Yrrl targets were mapped in the 
S96 yrrlA indicated yeast strains using ChlP-seq. (A) Heat map of overlapping peak regions across different alleles and conditions. Log 2 
fold changes [log2(FCs|] between normalized ChIP and control reads for each overlapping peak region were determined and compared 
across different yeast strains. A list of the overlapping peak regions and their log 2 (FC) values are provided in Supplemental Table 4. [B) 
Yrrl ChlP-seq signal tracks from three regions of the genome, including SNQ2 and YAR068w and between FMP52 and PMI40. The 
Y-axis is log 2 (FC), and the X-axis is the chromosome location with the chromosome and coordinates of each region noted. 
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Yrrl YJM789 in 4NQO, yrrl S96 177 SE in YPD, and yrrl S96 
I775E in 4NQO did bind known transcriptional targets,- for 
example, the promoter of SNQ2. Snq2 is an ABC plasma 
membrane transporter that is required for resistance to 
many chemicals, including 4NQO (Servos et al. 1993) and 
reactive oxygen species (ROS) (Ververidis et al. 2001), and 
is known to be regulated by Yrrl (Le Crom et al. 2002). 
For many loci, mutation of isoleucine 775 to glutamate in 
the S96 allele of Yrrl mimicked 4NQO treatment of 

Y rr i Y JW89 

Mapping Yrrl-binding regions across genetically 
diverged yeast strains 

The differences in phenotypes conferred by the YRR1 
alleles from the two strain backgrounds could be due to 
a number of possibilities. To determine how much the 
binding of Yrrl could change in both 4NQO and glycerol, 
we therefore mapped the binding sites of Yrrl in the five 
different strain backgrounds. Yrrl from AWRI1631, 
RM11, and YJM339 contains other sequence differences 
compared with S96 and YJM789 that alter binding to 
distinct gene targets. For each strain, targets were mapped 
in cells grown under three different conditions: rich 
medium (YPD), 4NQO, and glycerol. The sequencing 
reads were generated and mapped directly to their re- 
spective yeast genomes to avoid mapping biases (see the 
Materials and Methods,- Supplemental Table 5). Signifi- 
cant Yrrl peaks were scored with CisGenome (Ji et al. 
2008) using a variety of parameters (Supplemental Table 
6). For the results described below, we used a reasonably 
stringent threshold that provided a low false discovery 
rate (FDR,- 1 %),- qualitatively similar results were observed 
using a range of different thresholds (see the Supplemental 
Material for further discussion). 

Overall, we identified 125-411 binding regions in the 
different yeast strains and conditions (Table 1 ) ; in general, 



the highest number of total peaks was identified in cells 
grown in YPD [n = 41 1-256), and fewer were found in cells 
grown in the presence of 4NQO (n = 376-125) or glycerol 
[n = 387-130). Analysis of all of the genes associated 
with Yrrl -bound regions from all ChlP-seq experiments 
showed enrichment for genes involved in drug response 
even when cells were grown in the absence of drugs in 
rich medium or glycerol (Supplemental Tables 7, 8). 

Yrrl has different gene targets in different yeast 
backgrounds and conditions 

Analysis of Yrrl-binding regions revealed that there are 
extensive differences in both the different yeasts grown 
under the same condition and the same yeast grown 
under different conditions. To some extent, most peak 
regions bound by Yrrl changed between either strains or 
conditions. Examples of these variations are shown in 
Figure 6 A. 

The peak region in the intergenic space between YNG1 
and CYT1, which encode a subunit of the NuA3 histone 
acetyltransf erase and cytochrome cl, respectively, de- 
creased in AWRI1631, RM11, and YJM339 in 4NQO 
and glycerol media relative to S96 and YJM789. Cytl 
localizes to the mitochondria and has a role in heme 
biosynthesis. Varied binding patterns upstream of either 
gene could be partially responsible for the observed 
phenotype. An example is the Yrrl peak upstream of 
SNQ2, which is one of the strongest peaks in most strains 
and did not vary between different alleles of Yrrl seen in 
Figure 5B. The peak was strongest in AWRI1631 grown in 
YPD and YJM789 grown in glycerol, while peaks in other 
strains and conditions were about half as strong. The 
pattern of peaks upstream of GIS1 was similar across all 
conditions in AWRI1631, RM11, and YJM339 but varied 
greatly between S96 and YJM789. Gisl is a Zn finger TF 
with a jumonji demethylase domain regulating phospho- 



Table 1. Yrrl ChIP target and overlaps between strains 



Medium 


Strain 


Peaks 


YJM789 


S96 


AWRI1631 


RM11 


YJM339 


YPD 


YJM789 


277 


1 


6 


14 


5 


9 


YPD 


S96 


411 


21 




39 


7 


9 


YPD 


AWRI1631 


262 


11 


17 




11 


11 


YPD 


RM11 


387 


31 


9 


28 




31 


YPD 


YJM339 


259 


7 


17 


3 


11 




4NQO 


YJM789 


148 




33 


18 


9 


7 


4NQO 


S96 


125 


34 


10 


51 


44 


23 


4NQO 


AWRI1631 


138 


42 


33 




0 


4 


4NQO 


RM11 


131 


8 


40 


5 


11 


1 


4NQO 


YJM339 


376 


49 


64 


0 


0 


4 


Glycerol 


YJM789 


130 


59 


60 


89 


61 


61 


Glycerol 


S96 


144 


39 




79 


18 


14 


Glycerol 


AWRI1631 


138 


42 


0 




0 


0 


Glycerol 


RM11 


131 


47 


10 


29 




2 


Glycerol 


YJM339 


387 


49 


64 


0 


0 


4 



Number and percentage of unique peaks for each pairwise strain/condition combination. Off -diagonal entries are the number of peaks 
substantially different (log 2 normalized signal threshold) in the two-way comparison; e.g. in the cell at the intersection of row S96 (41 1) 
and column YJM339, 9% of the 411 binding regions in S96 are different in the two-way comparison of S96 and YJM339 under YPD 
growth condition. Along diagonal entries are the number of unique binding regions (shaded) in one strain but different from the other 
four strains under each condition; e.g. at the intersection of row S96 (411) and column S96, 2% of the 411 binding regions in S96 are 
different from other binding regions in the other four strains under YPD growth condition. 
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1 Kb Figure 6. Differential binding at distinct targets from 

five different strains in three different growth condi- 
tions. [A] Yrrl ChlP-seq signal tracks from YJM789 
(blue), S96 (red), AWRI1631 (green), RM11 (purple), and 
YJM339 (orange) were drawn for peak regions in YPD, 
4NQO, and glycerol media with the input tracks below 
each ChIP track. Genes with 5 ' ends downstream from 
peak regions ae indicated above the signal tracks, and 
corresponding regions on chromosomes are shown 
below. The black line represented the chromosome, 
and the scale is noted by a bar demonstrating 1 kb (X- 
axis). ORFs are boxed in black. The scale for each signal 
track for a condition was drawn to fit and is noted 
above the signal tracks (Y-axis). [B) Yrrl ChlP-seq signal 
tracks from YJM789 (blue) and S96 (red) from the three 
different conditions were drawn in the same scale for 

chr each peak region (Y-axis). The same schematic was used 

as described above to represent genes. Note that differ- 
ential Yrrl binding was observed for the different 
yeasts; in some cases, binding is similar under one 
condition but different in another (e.g., QCR2 in YPD 
compared with glycerol); in other cases, binding is 
different in all conditions (e.g., FLR1). 
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lipids levels (Oshiro et al. 2003), transcription of starva- 
tion genes, and DNA damage response (Jang et al. 1999). 
The GIS1 peak was similar across different strains in YPD 
but increased in YJM789 grown in glycerol and also in S96 
grown in 4NQO compared with other strains grown in 
the same conditions. The binding pattern between IPT1 
and SNF11 was similar to that observed for GIS1. Snfll 
is a component of the Swi/Snf chromatin remodeling 
complex (Treich et al. 1995) and Iptl, an inositolphospho- 
transf erase that has a role in drug resistance ( Halls trom 
and Moye-Rowley 2000). The binding pattern of the Yrrl 
upstream two genes encoding ribosomal proteins also 
was similar to the strongest binding observed in YJM339 
grown in 4NQO. Overall, the most variation in the number 
of peaks among the different yeast was observed in 4NQO 
(STD = 107.9) and glycerol (STD = 112.5) as compared with 
YPD (STD = 73.6). 



It is possible that the variation in Yrrl binding across 
all strains was due to an overall reduced level of Yrrl's 
ability to bind DNA in different strains. However, we do 
not believe that this is the case because we observed some 
peaks with similar levels among the different strains (e.g., 
the intergenic peak between YNG1 and CYT1 varied less 
than twofold for strains grown in YPD) and because we 
observed some very strong peaks in any given strain 
and condition (e.g., between IME1 and RPL43b). We also 
performed a systematic analysis, as described in the next 
section. 

Varied patterns of binding across different strains 
and conditions revealed by systematic analysis 

To compare all peaks among the different strains, we used 
the stringent criteria described in the Materials and 
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Methods and a pipeline similar to that described pre- 
viously (Kasowski et al. 2010; Zheng et al. 2010). We first 
identified the regions that bound Yrrl in the different 
strains and determined the normalized signals in these 
regions for each of these strains. We then analyzed the 
differences in signals as well as the percentage overlap of 
the binding regions. We first determined the correlation 
of signals in all Yrrl -bound regions between (1) the same 
strain grown in different conditions (Supplemental Fig. 8), 
(2) different strains in the same condition, and (3) as a 
control, the independent biological replicates from the 
same strain. The individual comparisons are presented 
in Supplemental Figures S4 and S5, and the results are 
summarized in Figure 7 A. For the Yrrl -bound regions, the 
correlations in normalized signals between different 
strains and different conditions were significantly lower 
than the correlations between biological replicates ( Wilcoxon 
rank-sum test, P = 0.02635 and P = 0.02042, respectively) 
(Fig. 6 A), indicating significant variation in binding among 
the different strains grown in the same condition and the 
same strain grown under different conditions. Overall, 
there was just as much variation in Yrrl targets between 
strains as there was between different conditions. 

For each binding target location, we also examined 
the differences in normalized binding signals for all five 
strains and all three conditions using the clustering 
diagram depicted in Figure 6B. As is evident in the figure, 
extensive Yrrl target and signal differences exist between 
each of the different strains irrespective of the growth 
conditions. However, some strains do appear more sim- 
ilar to one another under the same conditions (e.g., 
AWRI1631 and RM11 in 4NQO and RM11 and S96 in 
glycerol) (Supplemental Fig. S5). Interestingly, in at least 
two cases, the patterns appeared similar for different 
yeast exposed to different conditions (e.g., YJM789 and 
S96 grown in glycerol and 4NQO, respectively). These 
results are also evident from the pairwise correlation 
analyses (Supplemental Fig. S6) and suggest that the yeast 
TFs can switch their binding programs. 

We also determined the number of binding sites that 
were unique or substantially different in signal (more 
than twofold apart in normalized signal) between the 
different yeast strains grown under both the same and 
different conditions. Consistent with the heat map re- 
sults (Table 1), we observed extensive divergence in 
similar binding sites among the different yeast grown in 
the same condition as well as for the same yeast grown 
in different conditions. For example, YJM789 and S96 had 
33% and 40% unique/different binding sites (more than 
twofold) for cells grown on glycerol or in the presence of 
4NQO, respectively. To ensure that these differences 
were not due to threshold effects, two analyses were 
performed. First, we changed the threshold of fold differ- 
ences and found similar results (data not shown). Second, 
we examined the signals of peaks that were not found to 
be overlapping; potentially, many of these might lie just 
below the threshold. As shown in Supplemental Figure 
S7, although some peaks may be similar and fall below 
the threshold in one of the strains, a large number of 
peaks were very different, and the Yrrl ChIP signal was 
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Figure 7. Different binding targets of Yrrl among five yeast 
strains. Yrrl targets were mapped in five indicated yeast strains 
using ChlP-seq. {A) A summary of the pairwise correlations 
between peak signals (log read counts) in ChIP samples of two 
biological replicates of the same strain {left), the same strain 
under different conditions {middle), and different strains under 
the same condition {right). Correlations between different 
strains grown under the same conditions or the same strains 
grown under different conditions were significantly different 
from that observed for biological replicates (significance differ- 
ences between replicates and conditions [Wilcoxon rank-sum 
test, P = 0.026] or different strains [Wilcoxon rank-sum test, P = 
0.020]), whereas correlations between conditions and between 
strains were not significantly different (Wilcoxon rank-sum test, 
P = 0.4526). {B) Heat map of overlapping peak regions across 
different strains and conditions. Log 2 fold changes between 
normalized ChIP and control reads for each overlapping peak 
region were determined and compared across different yeast. 
A list of the overlapping peak regions and their log2(FC) values 
is provided in Supplemental Table 7. 

near background levels in one strain or the other. Overall, 
these results indicate that both the number of Yrrl 
targets and their binding signals often vary among differ- 
ent yeasts. 

Extensive variations in Yrrl localization were evident 
at different loci in the different yeast and conditions. 
Common overall patterns were evident based on the 
dendrogram above the heat map (Fig. 7B). In general, 
Yrrl from cells grown in glycerol had more similar 
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patterns of binding than that from cells grown in YPD or 
4NQO. The two Yrrl -binding patterns that were the 
most similar were in S96 grown in 4NQO and in 
YJM789 grown in glycerol. PDR3, SNQ2, PDR15, 
YOR1, PDR16, and PDR5, which were strongly associ- 
ated with Yrrl in all strains and conditions, were involved 
in drug response including Yrrl itself (Supplemental 
Table 8). 

From the peaks identified from the ChlP-seq analysis, 
we also performed a de novo motif search for Yrrl- 
binding sites (see the Materials and Methods) in all strains 
and growth conditions using MEME and recovered the 
consensus binding sites (Supplemental Fig. S8). These 
searches revealed distinct but related consensus binding 
sites among one another, which were also related to 
the published consensus (RCGGRDWTWD) (Zhu et al. 
2009). In particular, the consensus binding site from S96 
in 4NQO was very similar to the binding sites seen in 
S96, RM11, and YJM339 grown in glycerol. AWRI1631 
binding sites from glycerol showed some similarity to 
YJM789 grown in YPD. In many cases, when there was 
a similar deduced motif, there was also a similar Yrrl- 
binding pattern (e.g., AWRI1631 in 4NQO and glycerol 
media). However, overall, we did observe condition- and 
strain-specific variations, suggesting that the binding pref- 
erence may differ across different growth conditions. 



Discussion 

We identified and characterized a single protein (Yrrl) 
that affected growth in three different conditions and is 
polymorphic across diverse yeast strains. The polymor- 
phisms in YRR1 were outside the DNA-binding region, 
with the nonconservative substitutions located at the C 
terminus. One of the residues, at position 775, was likely 
the major reason for these different phenotypes, as 
changing it to a glutamate causes 4NQO resistance in 
an otherwise sensitive strain. Since conversion to an 
alanine does not eliminate 4NQO resistance in YJM789, 
we propose that other residues are responsible for the 
phenotypic variations as well. Without glutamate at 673, 
the Yrr 1 YIM789 could not confer 4NQO resistance without 
mutation of 775 to a glutamate. Unlike the I775E sub- 
stitution in Yrr I s96 , mutation of 673 to a glutamate could 
not convert the Yrr I s96 allele to a 4NQO-resistant allele. 
We further propose that these Yrrl amino acid differences 
lead to variations in association with protein partners/ 
cofactors such as kinases, which in turn results in 
selection of different binding targets. Since previous 
studies have shown in yeast and humans cells that TFs 
often have different partners at distinct loci that are 
important for binding (Kasowski et al. 2010; Zheng 
et al. 2010; Karczewski et al. 2011), it is likely that 
Yrrl has distinct partners at different targets, varying in 
the different strains and under the different conditions. A 
likely candidate is a serine/threonine kinase that phos- 
phorylates polymorphic sites in Yrrl in response to 
4NQO. The interaction of regulatory enzymes and their 
substrates is difficult to characterize and often requires 



overexpression or in vitro purification to monitor phys- 
ical interactions (Ruiz et al. 2011). It is interesting that 
Yrrl targets are sometimes more similar in different 
strains grown under different conditions then in either 
the same strain grown under different conditions or differ- 
ent strains grown in the same conditions. The simplest 
explanation for this observation is that different factors 
work with Yrrl under different strains and conditions and 
that some of these combinations are shared in those 
cases. It is likely the amino acid variations that are re- 
sponsible for these interaction differences, which is 
consistent with our demonstrating that coding sequence 
changes are primarily responsible for the phenotypic 
differences that we observed. 

Yrrl is polymorphic among different yeast strains, and 
we propose that this TF belongs to a class of regulatory 
proteins that we call "master variators." We suggest that 
master variators are proteins that have three characteris- 
tics: (1) They are polymorphic,- (2) they are regulatory 
proteins such as TFs, proteins kinases, and other modifi- 
cation enzymes,- and (3) they affect significant numbers of 
downstream targets. In contrast to deletions, study of 
natural variation would likely limit suppressor mutations 
that would compensate for severe phenotypes (Teng et al. 
2013). Variation in Yrrl and other master regulators can 
cause significant differences in phenotypes through alter- 
ation in regulatory programs rather than simply loss or 
reduction of function, as has been described previously 
(Chan et al. 2010). These alterations can improve one 
phenotype at a cost to another,- indeed, the Yrrl YJM789 
allele improves growth in the presence of 4NQO but 
reduces growth in glycerol even when the allele respon- 
sible for improved growth on glycerol is present. Varia- 
tion of Yrrl between ECU 18 and S288c affects the 
regulation of SNG1 during the fermentation process 
(Brion et al. 2013). We note that this concept is different 
from that reported by others in which loss of function 
causes a distinct phenotype. For Yrrl, we observed a 
phenotype transformation due to a single gene through 
molecular reprogramming. 

The mechanism by which Yrrl causes distinct pheno- 
types is not known in detail but may involve differences 
in many different binding targets. Yrrl binds upstream of 
many genes affecting mitochondrial function and drug 
response, and these targets differ in the different strain 
backgrounds (Supplemental Table 9). For example, Yrrl is 
differentially bound upstream of ISU1, QCR2, and TBS1 
(genes involved in mitochondrial function in S96 and 
YJM789) in some conditions, and in other conditions, it 
binds upstream of PDR5, SNG1, and FLR1 (genes in- 
volved in drug response) (Fig. 6B). Yrrl activity correlates 
with SNG1 and FLR1 expression (Le Crom et al. 2002; 
Teixeira et al. 2008; Brion et al. 2013; Kodo et al. 2013), 
and thus Yrrl is likely a direct activator of gene expres- 
sion of these targets. Yrrl binding in YJM789 at these and 
other regions may facilitate 4NQO resistance. Interest- 
ingly, Yrrl is thought to be a repressor as well as an 
activator because its activity can affect genes in different 
ways (Cui et al. 1998; Zhang et al. 2001; Le Crom et al. 
2002,- Lucau-Danila et al. 2003). Thus, its variation in 
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binding at different gene targets can alter expression in 
different ways,- as such, its control over variation as 
a master variator likely affects gene expression in very 
different ways and thus leads to the great extent of pheno- 
typic diversity. 

In conclusion, we demonstrate that variation in an 
important TF is responsible for major phenotypic changes 
in yeast. Previous studies have highlighted the elevated 
levels of nonsynonymous changes in TFs among different 
yeast strains (Gu et al. 2005). Since many human TFs and 
other regulatory proteins are also polymorphic, it is likely 
that some of these potential master variators may also be 
responsible for phenotypic differences, and thus this may 
be a general mechanism for mediating extensive regula- 
tory (e.g., expression and binding QTLs) and phenotypic 
variation in all eukaryotes. 

Materials and methods 

Yeast strains and plasmids 

All yeast strains used in this study are haploids and are described 
in Supplemental Table 10. S96 (isogenic to S288c) is a laboratory 
strain, YJM789 was derived from a clinical isolate from the 
lungs of an AIDS patient (McCusker et al. 1994; Wei et al. 2007), 
AWRI1631 was derived from a wine strain (Borneman et al. 
2008), and RMll-lb was derived from the vineyard diploid by 
removing HO and was sporulated, and the KAN R was removed as 
described (Gueldener et al. 2002). YJM339-2d was derived from 
a clinical isolate from a human bile tube (McCusker et al. 1994), 
and a stable haploid was generated as RMll-lb. pGS35 (KAN R ) 
and pGS36 (HYG R ) are centromeric plasmids that were a gift 
from Gavin Sherlock. All oligonucleotides used are described in 
Supplemental Table 11. YRR1 with the endogenous promoter 
and terminator were cloned by homologous recombination into 
the Xbal and Spel sites into pGS35. The entire transcriptional 
unit of YRR1 was selected, and the promoter contained four 
polymorphs, while downstream polymorphisms occured down- 
stream from the transcriptional stop of the YRR1 mRNA (Xu 
et al. 2009). Yeast were grown under standard conditions and 
transformed as outlined in Gietz and Schiestl (2007). YRR1 was 
tagged both at the C terminus with the 13xMyc tag marked with 
KAN R (Longtine et al. 1998) and at the N terminus by PCR- 
stitching the 13xMyc tag to the coding region of YRR1 YJM789 and 
YRR1 S96 and cloning into pSG35. Endogenous YRR1 was deleted 
via homologous recombination with HYG R (Goldstein and 
McCusker 1999), and plasmids were transformed into strains 
using standard lithium acetate transformation. Plasmids were 
maintained by addition of 200 |xg/mL G418 for selection. Site- 
directed mutagenic PCR was carried out with overlapping primers 
containing each of the mutations and flanking primers used for 
cloning. The 5' and 3' PCR products were then used as templates 
to PCR-stitch yrrl with the mutation, yrrl alleles were cloned into 
pGS35, and mutations were verified by sequencing the plasmid. 
Growth rates were calculated by measuring the OD 60 o of a yeast 
culture grown to log phase in YPD, washed, and suspended in 
30 mL of 4NQO or YPglycerol. Readings while the culture was in 
exponential growth were fitted to a log scale. The doubling time 
in hours was derived by the equation fitted to the growth curve. 

Immunoprecipitation 

Cells were grown in YPD to mid-log phase and then shifted to 
other medium. Next, 0.25 |xg/mL 4NQO dissolved in DMSO 
( Sigma- Aldrich) was added for 2.5 h to YPD, and/or cells were 



washed with YP and shifted to YP + glycerol for 8.5 h. Cells were 
lysed in the Fastprep for 20 sec at speed 5 and then ice for 5 min 
three times in 50 mM Tris HCl (pH 7.8), 150 mM NaCl, 4% 
glycerol, 1.5 mM Mg acetate, and 0.15% NP-40. Clarified extract 
(1.5 mg) was incubated with 35 julL of 50% anti-Myc EZ-view 
affinity gel (Sigma- Aldrich) for 2 h to overnight at 4°C. Washed 
immunoprecipitates were boiled with Lammeli sample buffer, 
separated on 4%-12% polyacrylamide gel, and transferred to 
nitrocellulose. Westerns were blotted with mouse anti-Myc 
9E10 (Sigma). To determine whether there was any background 
in the immunoprecipitations, a mock immunoprecipitation was 
carried out with Yrrl -Myc- tagged lysate and protein A-agarose 
beads. Lysate was treated with 15 julL of CIP overnight at 4°C, and 
then Yrrl -Myc was immunoprecipitated. 

4NQO response QTL mapping 

A total of 124 S96/YJM789 segregants (Mancera et al. 2008) was 
manually scored for growth on 4NQO medium on a scale of 0-5. 
Their specific genotypes at 55958 SNP marker loci were filtered 
(must have more than four strains of either parental genotype) 
and reduced to 5523 nonredundant genotype blocks. One- 
dimensional QTL scanning was performed using the R package 
qtl with EM method. The 95% Bayes credible interval was 
calculated for the highest peak on chromosome 15. The signif- 
icance threshold was above LOD 3.84 after 1000 permutations, 
and the peak above YRR1 was the only QTL that passed the 
threshold. The peak was represented by 10 A LOD rescaled to make 
the area under the curve equal 1 . 

ChlP-seq 

Cells were grown in YPD to mid-log phase, diluted to early log 
phase, and split into three cultures. YPD cultures were grown for 
2.5 h. 4NQO (0.25 |xg/mL) dissolved in DMSO (Sigma) was added 
for 2.5 h. The third culture was washed with YP and shifted to YP 
+ glycerol for 8.5 h. Thirty-five OD units were cross-linked with 
1% formaldehyde for 45 min before being washed once with 
water and flash-frozen. ChIP was carried out as described 
(Lefrancois et al. 2009) with the following modifications: The 
addition of glycine was omitted, and 35 julL of 50% anti-Myc EZ- 
view affinity gel was added (Sigma- Aldrich). Two-hundred nano- 
grams of chromatin was used as input controls. The construction 
of multiplexed barcoded libraries permitted the simultaneous 
sequencing of 24 samples in the Illumina HiSeq sequencers 
(Lefrancois et al. 2009). Libraries were sequenced on Illumina 
HiSeq as 100-base-pair (bp) paired-end reads. For different allelic 
ChIP of pMyc-Yrrl, there were three biological replicates for 
ChIP of each treatment and one replicate for input. The average 
generated reads were 3,357,478, and the alignment rate was 
99.17% for ChIP samples,- the average generated reads were 
3,134,818, and the alignment rate was 98.70% for input. ChIP 
from different strains with Yrrl -Myc tag was carried out as 
described (Aparicio et al. 1991 ). Two to three biological replicates 
of genetically divergent strain ChIP samples were sequenced on 
Illumina GAII as 32-bp single-end reads, generating an average of 
2.66 million post-filtered reads, with 64% mapping back to the 
genome; the input samples generated an average of 3.55 million 
post-filtered reads, with 77% mapping back to the genome. Nor- 
malized peak signals on log 2 scale for all strains and conditions were 
clustered by Euclidian distance using heatmap.2 function in R. 
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